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Summary. We compute numerically eigenvalues and eigenfunctions of the Lapla- 
cian in a three-dimensional hyperbolic space. Applying the results to cosmology, we 
demonstrate that the methods learned in quantum chaos can be used in other fields 
of research. 



1 Introduction 




V 




Fig. 1. Applications of the Laplacian. A drum (top left), electromagnetic waves 
(top middle), semiconductors resp. quantum mechanics (top right), and the universe 
(bottom). 



The Laplacian A, a second order differential operator, is of fundamental 
interest in several fields of mathematics and physics. One of the oldest exam- 
ples is a drum T>, see figure 1 top left. If one hits it, its membrane oscillates 
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and gives a sound. The vibrations of the drum are solutions of the Helmholtz 
equation (c = 1), 

(A-^)u(x,t)=0 VxeV, (1) 

where u is the displacement of the membrane. Fixing the membrane to its 
frame gives rise to Dirichlet boundary conditions, 

u(x,t)=0 yxedV. (2) 

A Fourier transformation, 

/oo 
v{x,u)e- iuJt du>, (3) 
-oo 

allows us to eliminate the time-dependence. The sound of the drum is deter- 
mined by the eigenvalues of the time-independent Helmholtz equation, 

(A + lu 2 )v{x,uj) = VieD, (4) 

which is nothing else than the eigenvalue equation of the (negative) Laplacian 
inside the domain T>. 

Another example is the propagation of electromagnetic waves in a two- 
dimensional system T>, see figure 1 top middle. Each component of the electric 
and magnetic field is again subject to the Helmholtz equation (1). Placing 
(super)conducting materials around V yields Dirichlet boundary conditions 
(2). Electromagnetic waves inside V can only be transmitted and received, if 
their frequencies arc in the spectrum of the Laplacian. 

Of special interest to this proceedings is quantum chaos which yields our 
next example, see figure 1 top right. A non-relativistic point particle moving 
freely in a manifold resp. orbifold V is described by the Schrodinger equation 

d H 2 

ih—<P(x, t) = A&(x, t) \fx e T>, (5) 

dt ' ' 2m y 1 w 

with appropriate boundary conditions on &D. Scaling the units to K = 2m = 1, 
and making the ansatz 

&(x,t) = ip(x)e-i Et , (6) 

gives the time-independent Schrodinger equation, 

{A + E)i>{x) = Vx e V. (7) 

The statistical properties of its spectrum and its eigenfunctions are a central 
subject of study in quantum chaos. 

In sections 8 and 9 we present some of the statistical properties of the 
solutions of the eigenvalue equation (7). But before, in sections 3-5 we in- 
troduce the three-dimensional hyperbolic system we are dealing with, and in 
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sections 6 and 7 we develop an efficient algorithm that allows us to compute 
the solutions numerically. 

From section 10 on we apply the results of the eigenvalue equation of the 
Laplacian to the universe, see figure 1 bottom, and compute the temperature 
fluctuations in the cosmic microwave background (CMB). This final example 
demonstrates that the methods learned in quantum chaos can be successfully 
used in other fields of research even though the physical interpretation can 
differ completely, e.g. metric perturbations and temperature fluctuations in 
cosmology instead of probability amplitudes in quantum mechanics. 

2 General statistical properties in quantum chaos 

Concerning the statistical properties of the eigenvalues and the eigenfunc- 
tions, we have to emphasise that they depend on the choice of the manifold 
resp. orbifold V. Depending on whether the corresponding classical system 
is integrable or not, there are some generally accepted conjectures about the 
nearest-neighbour spacing distributions of the eigenvalues in the semiclassical 
limit. The semiclassical limit is the limit of large eigenvalues, E — > oo. 

Unless otherwise stated, we use the following assumptions: The quantum 
mechanical system is desymmetrised with respect to all its unitary symme- 
tries, and whenever we examine the distribution of the eigenvalues we regard 
them on the scale of the mean level spacings. Moreover, it is generally believed 
that after desymmetrisation a generic quantum Hamiltonian corresponding to 
a classically strongly chaotic system possesses no degenerate eigenvalues. 

Conjecture 1 (Berry, Tabor [BT76]). If the corresponding classical system is 
integrable, the eigenvalues behave like independent random variables and the 
distribution of the nearest-neighbour spacings is in the semiclassical limit close 
to the Poisson distribution, i.e. there is no level repulsion. 

Conjecture 2 (Bohigas, Giannoni, Schmit [BGS84, BGS86]). If the corre- 
sponding classical system is chaotic, the eigenvalues are distributed like the 
eigenvalues of hermitian random matrices [Dys70, Meh91]. The corresponding 
ensembles depend only on the symmetries of the system: 

• For chaotic systems without time-reversal invariance the distribution of 
the eigenvalues approaches in the semiclassical limit the distribution of the 
Gaussian Unitary Ensemble (GUE) which is characterised by a quadratic 
level repulsion. 

• For chaotic systems with time-reversal invariance and integer spin the 
distribution of the eigenvalues approaches in the semiclassical limit the 
distribution of the Gaussian Orthogonal Ensemble (GOE) which is char- 
acterised by a linear level repulsion. 

• For chaotic systems with time-reversal invariance and half-integer spin the 
distribution of the eigenvalues approaches in the semiclassical limit the 
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distribution of the Gaussian Symplectic Ensemble (GSE) which is charac- 
terised by a quartic level repulsion. 

These conjectures are very well confirmed by numerical calculations, but sev- 
eral exceptions are known. Here are two examples: 

Exception 1 The harmonic oscillator is classically integrable, but its spec- 
trum is equidistant. 

Exception 2 The geodesic motion on surfaces with constant negative curva- 
ture provides a prime example for classical chaos. In some cases, however, 
the nearest-neighbour distribution of the eigenvalues of the Laplacian on these 
surfaces appears to be Poissonian. 

"A strange arithmetical structure of chaos" in the case of surfaces of constant 
negative curvature that are generated by arithmetic fundamental groups was 
discovered by Aurich and Steiner [AS88], see also Aurich, Bogomolny, and 
Steiner [ABS91]. (For the definition of an arithmetic group we refer the reader 
to [Bor69]). Deviations from the expected GOE-behaviour in the case of a 
particular arithmetic surface were numerically observed by Bohigas, Gian- 
noni, and Schmit [BGS86] and by Aurich and Steiner [AS89]. Computations 
showed [AS89, AS90], however, that the level statistics on 30 generic (i.e. non- 
arithmetic) surfaces were in nice agreement with the expected random-matrix 
theory prediction in accordance with conjecture 2. This has led Bogomolny, 
Georgeot, Giannoni, and Schmit [BGGS92], Bolte, Steil, and Steiner [BSS92], 
and Sarnak [Sar95] to introduce the concept of arithmetic quantum chaos. 

Conjecture 3 (Arithmetic Quantum Chaos). On surfaces of constant negative 
curvature that are generated by arithmetic fundamental groups, the distribu- 
tion of the eigenvalues of the quantum Hamiltonian approaches in the semi- 
classical limit the Poisson distribution. Due to level clustering small spacings 
occur comparably often. 

In order to carry out some specific numerical computations, we have to 
specify the manifold resp. orbifold V. We will choose V to be given by the 
quotient space T\H in the hyperbolic upper half space H where we choose T 
to be the Picard group, see section 5. 

3 The hyperbolic upper half-space 

Let 

H = {(x ,x 1 ,y) gR 3 ; y > 0} (8) 

be the upper half-space equipped with the hyperbolic metric of constant cur- 
vature — 1 
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dg2 = d _4±d4±dl 

y 2 

Due to the metric the Laplacian reads 

Q2 q2 q2 q 

A = y2 ^ + ^ + W ) - y d-y (10) 

and the volume element is 

dx dxidy 

dfi= - 3 . (11) 

The geodesies of a particle moving freely in the upper half-space are straight 
lines and semicircles perpendicular to the x -xi-plane, respectively, see figure 
2. 



y 




Fig. 2. Geodesies in the upper half-space of constant negative curvature. 

Expressing any point (x ,xi,y) e H as a Hamilton quaternion, z = x + 
ixi +]y, with the multiplication defined by i 2 = — 1, j 2 = — 1, ij +ji = 0, all 
motions in the upper half-space are given by linear fractional transformations 

z -fz — (az + b)(cz + d)~ ; a, b, c, G C, ad — be = I. (12) 

The group of these transformations is isomorphic to the group of matrices 

7=("J) eSL(2,C) (13) 
up to a common sign of the matrix entries, 

SL(2,C)/{±1} = PSL(2,C). (14) 

The motions provided by the elements of PSL(2, <D) exhaust all orientation 
preserving isometrics of the hyperbolic metric on H. 
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Remark 1. If one wants to avoid using quaternions, the point (x n ,x 1 ,y) E H 
can be expressed by (x, y) G C x E, with x = xo + ixi and y > 0. But then 
the linear fractional transformations look somewhat more complicated, 

,(ax + b)(cx + d) + acy 2 y , 

(x, y) ^ j{x, y) = ( : — [J , | — ^j. (15) 

v |ca; + dp + \cy\ z \cx + d\ l + \cy\ z ' 

In order to keep the notation simple, we mainly use quaternions. 



4 Topology 

The topology is given by the manifold rcsp. orbifold. An orbifold is a space that 
locally looks like a Euclidean space modulo the action of a discrete group, e.g. 
a rotation group. A manifold is an orbifold that locally resembles Euclidean 
space. 

A hyperbolic three-orbifold can be realised by a quotient 

r\H = {rz; z e H}, (16) 

where TL is the upper- half space and r is a discrete subgroup of the isometries 
on TL. The elements of r are 2 x 2-matrices whose determinants equal one. 
If the trace of an element is real, it is called hyperbolic, parabolic, or elliptic, 
depending on whether the absolute value of its trace is larger, equal, or smaller 
than two, respectively. If the trace of an element is not real, the element is 
called loxodromic. 

The action of r on TL identifies all T-equivalent points with each other. 
All the T-equivalent points of z give an orbit 

r z = { lz - 7 er}. (17) 

The set of all the orbits is the orbifold. If, except of the identity, _T contains 
only parabolic and hyperbolic elements, then r\H is a manifold. 

There exist infinitely many hyperbolic three-orbifolds. But opposed to the 
two-dimensional case there do not exist any deformations of hyperbolic three- 
orbifolds, because of the Mostow rigidity theorem [Mos73, Pra73]. 

Each hyperbolic three-orbifold has its specific volume which can be finite 
or infinite. The volumes are always bounded from below by a positive constant. 

If the volume of r\H is finite and if r does not contain any parabolic 
elements, then the orbifold is compact. If r does contain parabolic elements, 
then the orbifold has cusps and is non-compact. 

The hyperbolic three-manifold of smallest volume is unknown. Only a 
lower limit for the volume is proven to be [PrzOl] 



vol(TW) > 0.281. 



(18) 
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The smallest known hyperbolic three-manifold is the Weeks manifold [Wee85] 
whose volume is 

vo\(M) ~ 0.943. (19) 

The volume of the Weeks manifold is the smallest among the volumes of 
arithmetic hyperbolic three-manifolds [CFJR01]. 

In contrast, hyperbolic three-orbifolds are known which have smaller vol- 
umes, whereby the one with smallest volume is also unknown. A lower limit 
for the volume is [Mey88a, Mey88b] 

voi(r\w) > 8.2 • i(r 4 , (20) 

and the smallest known hyperbolic orbifold is the twofold extension of the 
tetrahedral Coxeter group CT(22). With 

r = CT(22)J = CT(22) 2 nIso + (W), (21) 

where Iso + (7i) are the orientation preserving isometries of 7i, the volume of 
the orbifold is 

vol(r\H) ~ 0.039. (22) 

The volume of this orbifold is the smallest among the volumes of arithmetic 
hyperbolic three-orbifolds [CF86] . 



5 The Picard group 

In the following we choose the hyperbolic three-orbifold r\H of constant 
negative curvature that is generated by the Picard group, 

r = PSL(2,Z[i]), (23) 

where Z[i] = % + iZ are the Gaussian integers. 

The Picard group is generated by the cosets of three elements, 



1 1\ 1 i\ 1 

oi r loir i o 



(24) 



which yield two translations and one inversion, 

2i— ► 2 + 1, zi-^z + i, z i ► z ^ . (25) 
The three motions generating P, together with the coset of the element 



i 

-i 



(26) 



that is isomorphic to the symmetry 

z = x + ]y izi = —x + jy, (27) 
can be used to construct the fundamental domain. 
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Definition 1. A fundamental domain of the discrete group r is a closed sub- 
set T C 7i with the following conditions: 

(i) T meets each orbit Fz at least once, 

(ii) if an orbit Fz does not meet the boundary of T it meets T at most once, 
(Hi) the boundary of T has Lebesgue measure zero. 

For the Picard group the fundamental domain of standard shape is 

F = i> = x Q + ix 1 + jy e H; -^<a; <i, < x 1 < i, \z\ > 1}, 

(28) 

with the absolute value of z being defined by \z\ = [x^ + x\ + y 2 )? , sec figure 
3. Identifying the faces of the fundamental domain according to the elements 




Fig. 3. The fundamental domain of the Picard group. 



of the group r leads to a realisation of the quotient space r\Tt, see figure 4. 

With the hyperbolic metric the quotient space r\TC inherits the structure 
of an orbifold that has one parabolic and four elliptic fix-points, 

1 . /3 1 .1 . [T .1 . [3 



JO °' * =J ' Z= 2 +J VI' Z= 2 +1 2 +J V2' Z = 1 2 +J VI- 

(29) 

The parabolic fix-point corresponds to a cusp at z = joo that is invariant 
under the parabolic elements 

and (;;). (30) 

Because of the hyperbolic metric the volume of the non-compact orbifold r\H 
is finite [Huml9], 
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*1 



*0 



Fig. 4. Identifying the faces of the fundamental domain according to the elements 
of the Picard group. 

vol(r\W) = = 0.30532186472 . . . , (31) 

where 

!/£Z[i]-{0} 

is the Dedckind zeta function. 



6 Maass waveforms 

We are interested in the smooth and square-intcgrablc eigenfunctions of the 
Laplacian in the orbifold r\H. A function being denned on the upper half- 
space that is invariant under all discrete linear fractional transformations, 

ip{jz) = ip{z) v 7 e r, (33) 

is called an automorphic function. Any automorphic function can be identified 
with a function living on the quotient space r\H, and vice versa any function 
being defined on the quotient space can be identified with an automorphic 
function living in the upper half-space. The functions we are interested in 
match the definition of Maass waveforms [Maa49a, Maa49b] . 

Definition 2. Let P be a discrete subgroup of the isometries T C PSL(2, C) 
containing parabolic elements. A function ip{z) is called a Maass waveform 
if it is an automorphic eigenfunction of the Laplacian that is smooth and 
square-integrable on the fundamental domain, 
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il>(z) € C°°(H), (34) 

ip{z) e L 2 (rW), (35) 
(zi + = VzeW, (36) 

^( 7 «) = ,p(z) v 7 e r, zeH. (37) 

Since the Maass waveforms are automorphic with respect to a discrete sub- 
group -T that contains parabolic elements, they are periodic in the directions 
perpendicular to the cusp. This allows to expand them into Fourier series. 
In case of the Picard group we have 

1>(z)=u{y)+ ]T a yK ik (2<K\p\y)e 2 ^ x , (38) 
/3ez[i]-{o} 

where $t@x is the real part of the complex scalar product (3x, and 

u(y)=h yl+ik + hyl - ik * k * > (39) 
[b 2 y + b 3 ylny if k = 0. 

Kik{t) is the K-Bessel function [Wat44] whose order is connected with the 
eigenvalue E by 

E = k 2 + 1. (40) 
If a Maass waveform vanishes in the cusp, 

lim ip(z) = 0, (41) 

z^joo 

it is called a Maass cusp form. 

According to the Roelcke-Selberg spectral resolution of the Laplacian 
[Sel56, Roe66, Hej83], its spectrum contains both a discrete and a contin- 
uous part. The discrete part is spanned by the constant eigenfunction ip ko 
and a countable number of Maass cusp forms ip kl , tpfa , ip k3 , . . . which we take 

to be ordered with increasing eigenvalues, = E ko < E kl < E k2 < £fc 3 < 

The continuous part of the spectrum E > 1 is spanned by the Eisenstein se- 
ries E(z, 1 + ik) which are known analytically [Kub73, EGM85]. The Fourier 
coefficients of the function Ak(1 + ik)E(z, 1 + ik) are given by 

b = A K (l + ik), h =A K (l-ik), a fj = 2 Y, I-!'"' ( 42 ) 

A,A»€Z[i] ^ 
A M =/3 

where /3 £ Z[i] - {0}, and 

A K (s) =4ir- s r( S )C K (s) (43) 



has an analytic continuation into the complex plane except for a pole at s = 1 . 
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Defining 



A K (l+ik)E(z,l + ik) 
y/Z\A K (l + ik)\ 



(44) 



the Eisenstein series vpf lscn (z) is real. 

Normalising the Maass cusp forms according to 



(45) 



we can expand any square integrable function 4> e L 2 (r\H) in terms of the 
Maass cusp forms and the Eisenstein series, [EGM98], 



is the Petersson scalar product. 

The discrete eigenvalues and their associated Maass cusp forms are not 
known analytically. Thus, one has to calculate them numerically. Previous 
calculations of eigenvalues for the Picard group can be found in [SG91, Hun96, 
GH96, Ste99]. By making use of the Hecke operators [Sta84, SG91, Hei92, 
HA93] and the multiplicative relations among the coefficients, Steil [Ste99] 
obtained a non-linear system of equations which allowed him to compute 
2545 consecutive eigenvalues. Another way is to extend Hejhal's algorithm 
[Hcj99] to three dimensions [The03]. Improving the procedure of finding the 
eigenvalues [The02], we computed 13950 consecutive eigenvalues and their 
corresponding eigenfunctions. 

7 Hejhal's algorithm 

Hejhal found a linear stable algorithm for computing Maass waveforms to- 
gether with their eigenvalues which he used for groups acting on the two- 
dimensional hyperbolic plane [Hej99], see also [SS02, Ave03] for some appli- 
cations. We extend this algorithm which is based on the Fourier expansion 
and the automorphy condition and apply it to the Picard group acting on the 
three-dimensional hyperbolic space. For the Picard group no small eigenvalues 
0<_E = fc 2 + l<l exist [Str94]. Therefore, k is real and the term u(y) in 
the Fourier expansion of Maass cusp forms vanishes. Due to the exponential 
decay of the K-Bessel function for large arguments, 




(46) 



where 




(47) 




for t 



oo 



(48) 
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and the polynomial bound of the coefficients [Maa49b] , 

a = O(\p\), |/?|-oo, (49) 
the absolutely convergent Fourier expansion can be truncated, 

1>{z)= Yl a p yK ik (2n\f3\y)e 2 ^ +[[e}}, (50) 

/3GZ[i]-{0} 
\0\<M 

if we bound y from below. Given e > 0, k, and y, we determine the smallest 
M = M(e, k, y) such that the inequalities 

2nMy>k and K ik (2wMy) < emax(K ik (t)) (51) 
hold. Larger y allow smaller M. In all remainder terms, 

[[s}}= ^yK- lk (27T\f3\y)e 2 ^, (52) 

/3ez[i]-{o} 

\0\>M 

the K-Bessel function decays exponentially in \/3\, and already the K-Bessel 
function of the first summand of the remainder terms is smaller than e times 
most of the K-Bessel functions in the sum of (50). Thus, the error [[e]] does 
at most marginally exceed e. The reason why [[e]] can exceed e somewhat is 
due to the possibility that the summands in (50) cancel each other, or that 
the coefficients in the remainder terms are larger than in (50). By a finite 
two-dimensional Fourier transformation the Fourier expansion (50) is solved 
for its coefficients 

a^yK ik (27r\j\y) = — !— £ ${x + jy) e" 2 ^ +[[e}], (53) 
^ H > xexp] 

where 7 e Z[i] — {0}, and X[i] is a two-dimensional equally distributed set of 
(2Q) 2 numbers, 



h = -Q + h-Q + l---,Q-lQ-h J = 0,1}, 

(54) 



with 2Q > M+ | 7 |. 

By automorphy we have 

#*)=#0> (55) 
where z* is the T-pullback of the point z into the fundamental domain T, 

z * = 1Z , 7 e r, z* e J". (56) 
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apy*K ik (2ir\f3\y*)e- 



27riSR/3x* 



-M, (57) 



/3ez[i]-{o} 

\P\<M 



where y* is always larger or equal than the height yo of the lowest points of 
the fundamental domain T , 



Vo 



(58) 



allowing us to replace M(e, k, y) by M = M(e, k, yo). 

Choosing y smaller than y the T-pullback z z* of any point into the 
fundamental domain T makes at least once use of the inversion z i— > —z , 
possibly together with the translations zhz|1 and znz|i. This is called 
implicit automorphy, since it guarantees the invariance tp(z) = tp(— z^ 1 ). The 
conditions ip(z) = ip(z + 1) and i[)(z) — ip(z + i) are automatically satisfied 
because of the Fourier expansion. 

Making use of the implicit automorphy by replacing tp(x+jy) in (53) with 
the right-hand side of (57) gives 



a^yK ife (27r|7|y) 



(2Q) 2 



J2 E apy*K ik (27r\f3\y*)c 2 ^e- 



27T1?R.JX 



+ [[2e}}, (59) 



xGX[i] /3eZ[i]-{0} 
\0\<M O 



which is the central identity in the algorithm. 

The symmetry in the Picard group and the symmetries of the fundamental 
domain imply that the Maass waveforms fall into four symmetry classes [Stc99] 
named D, G, C, and H, satisfying 



D 

G 
C 
H 



ip(x + jy) = ip{\x + jy) = ip(—x + jy), 
ip(x+ ]y) = ip(ix + jy) = -ip(-x + jy), 
i>(x + jy) = -ip(ix + jy) = ip(-x + jy), 
ip(% + jy) = -ip(ix + jy) = -ip(-x + jy), 



(60) 
(61) 
(62) 
(63) 



respectively, see figure 5, from which the symmetry relations among the coef- 
ficients follow, 



D 

G 
C 
H 



ap 
ap 
ap 
ap 



a ifj = 

Oi/3 = 
-ciip 



ap, 
-ap, 
= ap, 
= ~ap- 



(64) 
(65) 
(66) 
(67) 
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XQ 





XQ 



XQ 



Fig. 5. The symmetries D, G, C, and H from top left to bottom right. 



Defining 



where s a p is given by 



c S ((3,x)= ]T S ^e Mte , 

ctGSh 



and 



a e S/3 



{0, i/3, -0, -i0, 0, i~0, -0, -10} if # {0, 10, -0, -10}, 



K {0,i0,- 0,-10} else, 
the Fourier expansion (38) of the Maass waveforms can be written 
ip{z) = u{y) + a (i yK {k {2ir\0\y)cs{0,x), 

~;[i]-{o> 



where the tilde operator on a set of numbers is defined such that 
Id, (J § x = X, and P| § x = 



holds. 
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Forgetting about the error [[2e]] the set of equations (59) can be written 

as 

V 1 p(k,y)ap = 0, 7 eZ[i]-{0}, (73) 

/3ez[i]-{o} 

\P\<M 

where the matrix V = (Vyp) is given by 

V 70 (k,y)=yK ik (2ir\i\y)6 70 --^ ]T y*K ik (2n\(3\y*) cs(/3, x*) cs(- 7 , x). 

x£X[i] 

(74) 

Since y < yo can always be chosen such that Kik(2n\j\y) is not too small, the 
diagonal terms in the matrix V do not vanish for large \ j\ and the matrix is 
well conditioned. 

We are now looking for the non-trivial solutions of (73) for 1 < \-f\ < M 
that simultaneously give the eigenvalues E = k 2 + 1 and the coefficients ap. 
Trivial solutions are avoided by setting the first non- vanishing coefficient equal 
to one, a a = 1, where a is 1, 2 + i, 1, and 1 + i, for the symmetry classes D, 
G, C, and H, respectively. 

Since the eigenvalues are unknown, we discretise the fc-axis and solve for 
each fc-value on this grid the inhomogeneous system of equations 

V 1 p(k,y* 1 )ap = -V ia (k,y* 1 ), 1 < | 7 | < M 0) (75) 

/3GZ[i]-{0,a} 
\/3\<M 

where < yo is chosen such that Kik(2n\'-f\y : ^ 1 ) is not too small for 1 < 
|7| < M . A good value to try for is given by 27rM y #1 = k. 
It is important to check whether 

Vi0(k,y* 2 )ap, 1<|7|<M , (76) 

0£Z[i]-{O} 
\P\<M 

vanishes where y# 2 is another y value independent of y* 1 . Only if all g 7 vanish 
simultaneously the solution of (75) is independent of y. In this case E = k 2 + l 
is an eigenvalue and the a^'s are the coefficients of the Fourier expansion of 
the corresponding Maass cusp form. 

The probability to find a k- value such that all g 1 vanish simultaneously is 
zero, because the discrete eigenvalues are of measure zero in the real numbers. 
Therefore, we make use of the intermediate value theorem where we look for 
simultaneous sign changes in g 7 when k is varied. Once we have found them in 
at least half of the g 7 's, we have found an interval which contains an eigenvalue 
with high probability. By some bisection and interpolation we can see if this 
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interval really contains an eigenvalue, and by nesting up the interval until its 
size tends to zero we obtain the eigenvalue. 

It is conjectured [BGGS92, BSS92, Bol93, Sar95] that the eigenvalues of 
the Laplacian to cusp forms of each particular symmetry class possess a spac- 
ing distribution close to that of a Poisson random process, see conjecture 3. 
One therefore expects that small spacings will occur rather often (due to level 
clustering). In order not to miss eigenvalues which lie close together, we have 
to make sure that at least one point of the fc-grid lies between any two suc- 
cessive eigenvalues. On the other hand, we do not want to waste CPU time if 
there are large spacings. Therefore, we use an adaptive algorithm which tries 
to predict the next best k-vahie of the grid. It is based on the observation that 
the coefficients ap of two Maass cusp forms of successive eigenvalues must dif- 
fer. Assume that two eigenvalues lie close together and that the coefficients of 
the two Maass cusp forms do not differ much. Numerically then both Maass 
cusp forms would tend to be similar - which contradicts the fact that different 
Maass cusp forms are orthogonal to each other with respect to the Petersson 
scalar product 

{ij>k i ,1>k i )=0, iiki^kj. (77) 

Maass cusp forms corresponding to different eigenvalues are orthogonal be- 
cause the Laplacian is an essentially self-adjoint operator. Thus, if successive 
eigenvalues lie close together, the coefficients ap must change fast when vary- 
ing k. In contrast, if successive eigenvalues are separated by large spacings, 
numerically it turns out that often the coefficients change only slowly upon 
varying k. Defining 

~ af3 = —= = ^===, 1<\P\<M Q , (78) 
/2^ 7 ez;[i]-{o} Kl 

V |7|<M 

our adaptive algorithm predicts the next fc-value of the grid such that the 
change in the coefficients is 

Mfcoid)-M fc ncw)| 2 - 0.04. (79) 

!3ez[i\-{o} 

\/3\<M 

For this prediction, the last step in the fc-grid together with the last change in 
the coefficients is used to extrapolate linearly the choice for the next fc-value 
of the grid. 

However the adaptive algorithm is not a rigorous one. Sometimes the pre- 
diction of the next fc-value fails so that it is too close or too far away from 
the previous one. A small number of small steps does not bother us unless 
the step size tends to zero. But, if the step size is too large, such that the 
left-hand side of (79) exceeds 0.16, we reduce the step size and try again with 
a smaller fc- value. 

Compared to earlier algorithms, our adaptive one tends to miss signifi- 
cantly less eigenvalues per run. 
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8 Eigenvalues 

We have computed 13950 consecutive eigenvalues and their corresponding 
eigenfunctions of the (negative) Laplacian for the Picard group. The smallest 
non-trivial eigenvalue is E = k 2 + 1 with k = 6.6221193402528 which is in 
agreement with the lower bound E > -|- [Str94] . Table 1 shows the first few 
eigenvalues of each symmetry class and table 2 shows some larger ones. The 
eigenvalues listed in table 1 agree with those of Steil [Ste99] up to five decimal 
places. 

One may ask whether we have found all eigenvalues. The answer can be 
given by comparing our results with Weyl's law. Consider the level counting 
function (taking all symmetry classes into account), 

N(k) = #{ i | h < k}, (80) 

(where the trivial eigenvalue, E = 0, is excluded), and split it into two parts 

N(k) = N(k) + N fluc (k). (81) 

Here TV is a smooth function describing the average increase in the number of 
levels, and Nfi uc describes the fluctuations around the mean such that 

1 f K 

lim — / N fluc (k)dk = 0. (82) 

The average increase in the number of levels is given by Weyl's law [Weyl2, 
Ava56] and higher order corrections have been calculated by Matthies [Mat95]. 
She obtained 

N(k) = ^^k 3 + a 2 k log k + a 3 k + a 4 + o(l) (83) 

with the constants 

«2 = (84) 
a 3 = log 2 + I \ogir - \o g r{\) + | log(2 + V3) + f], (85) 
a 4 = -§. (86) 

We compare our results for N(k) with (83) by defining 

N fluc (k) = N(k)~N(k). (87) 

Nfiuc fluctuates around zero or a negative integer whose absolute value gives 
the number of missing eigenvalues, see figure 6. Unfortunately, our algorithm 
does not find all eigenvalues in one single run. In the first run it finds about 
97% of the eigenvalues. Apart from very few exceptions the remaining eigen- 
values are found in the third run. To be more specific, we plotted Nfi uc de- 
creased by i, because N(k) — N(k) is approximately \ whenever E = k 2 + 1 
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Table 1. The first few eigenvalues of the negative Laplacian for the Picard group. 
Listed is k, related to the eigenvalues via E — k 2 + 1. 



D 


G 


c 


H 


8.55525104 




6.62211934 




11.10856737 




10.18079978 




12.86991062 




12.11527484 


12. 


14 07966049 




12 87936900 




15 34827764 




14 14833073 




15.89184204 




14.95244267 


14. 


17.33640443 




16.20759420 




17.45131992 


17.45131992 


16.99496892 


16. 


17 77664065 




17 86305643 


17, 


19.06739052 




18.24391070 




19.22290266 




18.83298996 




19.41119126 




19.43054310 


19. 


20 00754583 




20 30030720 


20. 


20 70798880 


20 70798880 


20 60686743 




20.81526852 




21.37966055 


21. 


21.42887079 




21.44245892 




22.12230276 




21.83248972 


21. 


22.63055256 




22 58475297 


22. 


22.96230105 


22.96230105 


22.85429195 




23.49617692 




23.49768305 


23. 


23.52784503 




23.84275866 




23 88978413 


23 88978413 


23 89515755 

— ■ 1 . ! ' ■ 1 1 1 ' 1 I I 


23. 


24.O4O01664 




Q/i /loioooon 
24.4zlo38zy 


z4. 


24.57501426 




25.03278076 


25. 


24.70045917 




25.42905483 




25.47067539 




25.77588591 


25. 


25.50724616 




26.03903968 




25.72392169 


25.72392169 


26.12361823 


26. 


25.91864376 


25.91864376 


26.39170209 




26.42695914 




27.07065195 


27. 


27.03326136 




27.16341524 


27. 


27.14291906 




27.26799477 


27. 


27.14498438 


27.14498438 


27.89811315 





12.11527484 



14.95244267 



21.37966055 



23.49768305 
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Table 2. Some consecutive large eigenvalues of the negative Laplacian for the Picard 
group. Listed is k, related to the eigenvalues via E — k 2 + 1. 



D 




G 




c 




H 




139 


.65419675 


139 


.65419675 


139 


.66399548 


139 


.66399548 


139 


.65434417 


139 


.65434417 


139 


.66785333 


139 


.66785333 


139 


.65783548 


139 


.65783548 


139 


.66922266 


139 


.66922266 


139 


661 04047 


139 


66104047 


139 


67870460 


139 


67870460 


139 


67694018 






139 


.68234200 


139 


.68234200 


139 


.68162707 


139 


.68162707 


139 


.68424704 


139 


.68424704 


139 


.68657976 






139 


.69369972 


139 


.69369972 


139 


71803029 


139 


71803029 


139 


69413379 


139 


69413379 


139 


72166907 


139 


72166906 


139 


69657741 


139 


69657741 


139 


.78322452 


139 


.78322452 


139 


.73723373 


139 


.73723373 


139 


.81928622 


139 


.81928622 


139 


.73828541 


139 


.73828541 


139 


.81985670 


139 


.81985670 


139 


.74467774 


139 


.74467774 


139 


.82826034 


139 


.82826034 


139 


75178180 


139 


75178180 


139 


.84250751 






139 


.75260292 


139 


.75260292 


139 


.87781072 


139 


.87781072 


139 


.79620628 


139 


.79620628 


139 


.87805540 






139 


.80138072 


139 


.80138072 


139 


88211647 


139 


88211647 


139 


.81243991 


139 


.81243991 


139 


91782003 


139 


91782003 


139 


.81312982 


139 


.81312982 


139 


.91893517 






139 


.82871870 


139 


.82871870 


139 


.92397167 


139 


.92397167 


139 


.86401372 


139 


.86401372 


139 


.92721861 


139 


.92721861 


139 


.86461581 


139 


.86461581 


139 


93117207 

• 1 ' J 1 1 1 — V 1 1 


139 


93117207 


139 


89407865 
i . ' i.i i "ii 


139 


89407865 


i on 
139 


no i a now 
.931492/ / 


i on 
139 


no i /t noTT 
.9^1492 ( 1 


i on 
139 


onm a 777 
.89914/ 1 i 


i on 
139 


onm a 777 
.89914 1 1 l 


139 


.94067283 






139 


.90090849 


139 


.90090849 


139 


.94396890 


139 


.94396890 


139 


.91635302 


139 


.91635302 


139 


.95074070 






139 


.94071729 


139 


.94071729 


139 


.95124805 


139 


.95124805 


139 


.95080198 


139 


.95080198 


139 


.99098324 


139 


.99098324 


139 


.97043676 


139 


.97043676 


140 


.00011792 


140 


.00011792 


140 


.00409202 


140 


.00409202 


140 


.00109753 


140 


.00109753 


140 


.02733151 


140 


.02733151 


140 


.00626902 


140 


.00626902 


140 


.04198905 


140 


.04198905 


140 


.00827516 


140 


.00827516 


140 


.04799273 


140 


.04799273 


140. 


.01679122 


140 


.01679122 


140 


.05764030 


140. 


.05764030 
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is an eigenvalue. A plot indicating that Nf [uc fluctuates around zero is shown 
in figure 7 where we plotted the integral 

1 f K 

I(K) = — J N fluc (k)dk. (88) 
Desymmetrising the spectrum yields Weyl's law to be 




5000 10000 15000 



Fig. 6. Nfi uc (ki) as a function of N(ki) = i fluctuating around zero. 



0.5 

i 



Fig. 7. J as a function of K showing that / g ^°°> 0. 



N(k)--^ L ^ + 0(k") (89) 

for each symmetry class. Looking at table 1 it seems somehow surprising 
that there arc not equally many eigenvalues listed for each symmetry class. 
Especially in the symmetry classes G and H there seem to be much less 
eigenvalues than in the symmetry classes D and C. Indeed, as was shown by 
Steil [Ste99] , there occur systematic degenerated eigenvalues between different 
symmetry classes. 

Theorem 1 (Steil [Ste99]). If E = k 2 + 1 is an eigenvalue corresponding 
to an eigenfunction of the symmetry class G resp. H ; then there exists an 
eigenfunction of the symmetry class D resp. C corresponding to the same 
eigenvalue. 
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Based on our numerical results we conjecture [The03]: 

Conjecture 4- Taking all four symmetry classes together, there are no degen- 
erate eigenvalues other than those explained by Steil's theorem. Furthermore, 
the degenerate eigenvalues which are explained by Steil's theorem occur only 
in pairs of two degenerate eigenvalues. They never occur in sets of three or 
more degenerate eigenvalues. 

Looking at the semiclassical limit, E — > oo, we finally find that almost all 
eigenvalues are two-fold degenerated, see e.g. table 2, which is an immediate 
consequence of Weyl's law, Steil's theorem, and conjecture 4. This means that 
as E — > oo 

^{non-degenerate eigenvalues < E} 
#{all eigenvalues < E} 

Finally, we remark that the distribution of the eigenvalues of each individ- 
ual symmetry class agrees numerically with conjecture 3, see [Ste99, The03]. 



9 Eigenfunctions 

Concerning the eigenfunctions of the Laplacian, it is believed that they behave 
like random waves. The conjecture of Berry [Bcr77] predicts for each eigen- 
function in the semiclassical limit, E — > oo, a Gaussian value distribution, 

1 _ u 2 

dp(u) = — c 2^" du, (91) 

V27T(7 



inside any compact regular subregion F of T . This means that 

hm — — r = 1 92) 

fa d P (u) 

holds with variance 

* 2 = ^/ F l^)! 2 ^ (93) 

for any — oo < a < b < oo, where X[a,b] ^ s the indicator function of the interval 
[a,b]. Figure 8 shows the value distribution of the 80148th 3 eigenfunction 
corresponding to the eigenvalue E = k 2 + 1 with k = 250.0018575195 inside 
a small subregion 

F = {z = x + iy; -0.43750 < x < -0.31435, 

0.06250 < X! < 0.18565, (94) 
1.10000 < y < 1.34881}. 

Our numerical data agree quite well with Berry's conjecture, providing nu- 
merical evidence that the conjecture holds. A plot of the eigenfunction inside 
the region F is given in figure 9. 

3 The number 80148 was determined approximately using Weyl's law (83). 
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-0.2 w 0.2 

Fig. 8. In the top left figure the value distribution of the eigenfunction corresponding 
to the eigenvalue E = k 2 + 1 with k — 250.0018575195 inside the region F is shown 
as ( + ). The solid line is the conjectured Gaussian. In the figure on the top right 
the dashed line is the integrated value distribution of the eigenfunction which is 
nearly indistinguishable from the integrated Gaussian (solid line). The figure on the 
bottom right displays a detailed magnification of the integrated value distribution 
showing that it lies slightly below the integrated Gaussian. 

10 An application to cosmology 

In the remaining sections we apply the eigenvalues and eigenfunctions of the 
Laplacian to a perturbed Robertson- Walker universe and compute the tem- 
perature fluctuations in the cosmic microwave background (CMB). 

The CMB is a relic from the primeval fireball of the early universe. It is 
the light that comes from the time when the universe was 379 000 years old 
[BHH+03]. It was predicted by Gamow in 1948 and explained in detail by 
Peebles [Pee65]. In 1978, Penzias and Wilson [PW65] won the Nobel Prize 
for Physics for first measuring the CMB at a wavelength of 7.35 cm. Within 
the resolution of their experiment they found the CMB to be completely 
isotropic over the whole sky. Later with the much better resolution of the 
NASA satellite mission Cosmic Background Explorer (COBE), Smoot et al. 
[SBK+92] found fluctuations in the CMB which are of amplitude 10~ 5 relative 
to the mean background temperature of To = 2.725 K, except for the large 
dipole moment, see figure 10. These small fluctuations serve as a fingerprint 
of the early universe, since the temperature fluctuations are related to the 
density fluctuations at the time of last scattering. They show how isotropic 
the universe was at early times. In the inflationary scenario the fluctuations 
originate from quantum fluctuations which are inflated to macroscopic scales. 
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Fig. 9. A plot of the eigenfunction corresponding to the eigenvalue E — k 2 + 1 with 
k = 250.0018575195 inside the region F. 



Due to gravitational instabilities the fluctuations grow steadily and give rise 
to the formation of stars and galaxies. 

The theoretical framework in which the CMB and its fluctuations are 
explained is Einstein's general theory of relativity [Einl5a, Einl5b, Einl5c, 
Einl6, Einl7]. Thereby a homogeneous and isotropic background given by 
a Robertson- Walker universe [Fri22, Fri24, Lem27] is perturbed. The time- 
evolution of the perturbations can be computed in the framework of linear 
perturbation theory [Lif46, Bar80]. 

An explanation for the presence of the CMB is the following, see also fig- 
ure 11: We live in an expanding universe [Hub29]. At early enough times the 
universe was so hot and dense that it was filled with a hot plasma consisting of 
ionised atoms, unbounded electrons, and photons. Due to Thomson scattering 
of photons with electrons, the hot plasma was in thermal equilibrium and the 
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Fig. 10. Sky maps of the temperature fluctuations in the CMB as observed by the 
NASA satellite mission COBE. The sky map on the top shows the dipole anisotropy 
after the mean background temperature of To = 2.725 K has been subtracted. The 
amplitude of the dipole anisotropy is about 3mK. Also subtracting the dipole yields 
the sky map in the middle. One sees the small temperature fluctuations whose 
amplitude is roughly 30 /xK. But one also sees a lot of foreground contamination 
along the equator that comes from nearby stars in our galaxy. After removing the 
foregrounds one finally gets the sky map on the bottom showing the temperature 
fluctuations in the CMB. Downloaded from [GH03]. 
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13.8 Gy 



379 000 y 



today 



last scattering 



2.725 K 



neutral atoms 
& 

free photons 



3000 K 




plasma p+ e ~ y 



thermal equilibrium 



T T 



Fig. 11. The expanding universe. At the time of last scattering occured a phase 
transition from an opaque to a transparent universe. 




mean free path of the photons was small, hence the universe was opaque. Due 
to its expansion the universe cooled down and became less dense. When the 
universe was around 379 000 years old, its temperature T has dropped down 
to approximately 3000 K. At this time, called the time of last scattering, the 
electrons got bound to the nuclei forming a gas of neutral atoms, mainly hy- 
drogen and helium, and the universe became transparent. Since this time the 
photons travel freely on their geodesies through the universe. At the time of 
last scattering the photons had an energy distribution according to a Planck 
spectrum with temperature of nearly 3000 K. The further expansion of the 
universe redshifted the photons such that they nowadays have an energy dis- 
tribution according to a Planck spectrum with temperature of T = 2.725 K. 
This is what we observe as the CMB. 

Due to the thermal equilibrium before the time of last scattering the CMB 
is nearly perfectly isotropic, but small density fluctuations lead to small tem- 
perature fluctuations. The reason for the small temperature fluctuations comes 
from a variety of effects. The most dominant effects are the gravitational red- 
shift that is larger in the directions of overdense regions, the intrinsic temper- 
ature fluctuations, and the Doppler effect due to the velocity of the plasma. 



11 Robertson- Walker universes 



Assuming a universe whose spatial part is locally homogeneous and isotropic, 
its metric is given by the Robertson- Walker metric, 



ds 2 = dt 2 - a 2 (t)fijdx % dx j , 



(95) 
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where we use the Einstein summation convention. Notice that we have changed 
the notation slightly. Instead of the quaternion z for the spatial variables, we 
now write x = x a + ixi + ix 2 . "fij is the metric of a homogeneous and isotropic 
three-dimensional space, and the units are rescaled such that the speed of 
light is c = 1. Introducing the conformal time dr\ — -M* we have 

ds 2 = a 2 (r/) [drj 2 - j^dx'dx 1 ] , (96) 

where a(t]) = a(t(rj)) is the cosmic scale factor. 

With the Robertson- Walker metric the Einstein equations simplify to the 
Friedmann equations [Fri22, Fri24, Lem27]. One of the two Friedmann equa- 
tions reads 

a» + Ka 2 = ^T$a* + ±Aa A , (97) 

and the other Friedmann equation is equivalent to local energy conservation. 
a' is the derivative of the cosmic scale factor with respect to the conformal 
time rj. k is the curvature parameter which we choose to be negative, k = — 1. 
G is Newton's gravitational constant, T{f is the energy-momentum tensor, and 
A is the cosmological constant. 

Assuming the energy and matter in the universe to be a perfect fluid 
consisting of radiation, non-relativistic matter, and a cosmological constant, 
the time-time component of the energy-momentum tensor reads 

Tg = El { V )+e m {ri), (98) 
where the energy densities of radiation and matter scale like 

sM=eM(^) 4 and £m( , H£mM (^) 3 . (99) 

Here 770 denotes the conformal time at the present epoch. 

Specifying the initial conditions (Big Bang!) o(0) = 0, a'(0) > 0, the 
Friedmann equation (97) can be solved analytically [AS01], 

where V{r\) denotes the Weierstrass P-function which can numerically be eval- 
uated very efficiently, see [AS64] , by 

1 00 

Hri) = -P(V,92,93) = + V c nV 2n - 2 (101) 



V 



2 1 A^i 

n=1 



with 
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„ n-2 



92 



93 



for n > 4. 
(102) 



The so-called invariants 32 and 33 are determined by the cosmological param- 
eters, 



92 



Sl A Sl r 1 



m 12 



To' 93 — T 



1 o A n T 1 (Ia^L 1 



with 



f2 T 

f*c = 

Q A = 



6 ST* 16 STi 216' 



87rGe r (r? ) 



7777' ( 103 ) 



87rG , e m (?7o) 

3# 2 M ' 
1 

ff 2 Ma 2 (%)' 
yl 



where 



3iJ 2 (77o)' 



a 2 (,7) 



(104) 
(105) 
(106) 
(107) 

(108) 



is the Hubble parameter. 

Because of the homogeneity and isotropy, nowhere in the equations of a 
Robertson- Walker universe appears the Laplacian. 



12 Perturbed Robertson- Walker universes 

The idealisation to an exact homogeneous and isotropic universe was essential 
to derive the spacetime of the Robertson- Walker universe. But obviously, we 
do not live in a universe which is perfectly homogeneous and isotropic. We 
see individual stars, galaxies, and in between large empty space. Knowing the 
spacetime of the Robertson- Walker universe, we can study small perturbations 
around the homogeneous and isotropic background. Since the amplitude of the 
large scale fluctuations in the universe is of relative size 10~ 5 [SBK + 92], we 
can use linear perturbation theory. In longitudinal gauge the most general 
scalar perturbation of the Robertson- Walker metric reads 

ds 2 = a 2 (?7)[(l + 2<P)dr, 2 - (1 - 2<P)- Hj dx i dx j ], (109) 

where <P = <l>(r], x) and \P = ^(77, x) are functions of spacetime. 
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Assuming that the energy and matter density in the universe can be de- 
scribed by a perfect fluid, consisting of radiation, non-relativistic matter, and 
a cosmological constant, and neglecting possible entropy perturbations, the 
Einstein equations reduce in first order perturbation theory [MFB92] to 

& = (110) 
+ 3H(1 + c 2 )<P' - c 2 A<P + (2H 1 + (1 + 3c 2 )(H 2 + 1))<2> = 0, (111) 

where H = ^ and c 2 = (3 + f^-) -1 are given by the solution of the 
non-perturbed Robertson- Walker universe. In the partial differential equation 
(111) the Laplacian occurs. If the initial and the boundary conditions of <P are 
specified, the time-evolution of the metric perturbations can be computed. 
With the separation ansatz 

$(ri,x) =Y J fk{ri)Mx) + J dkfkWM*), ( 112 ) 

k 

where the ip k are the eigenfunctions of the negative Laplacian, and the E k are 
the corresponding eigenvalues, 

-Ail) k (x) =E k Mx), (H3) 

(111) simplifies to 

fib) + 3H(1 + cl)f k {r,) + {c 2 s E k + 2H> + (1 + 3c 2 s )(H 2 + l))/ fc („) = 0. 

(114) 

The ODE (114) can be computed numerically in a straightforward way, and 
we finally obtain the metric of the whole universe. This gives the input to the 
Sachs- Wolfe formula [SW67] which connects the metric perturbations with 
the temperature fluctuations, 

ST 3 f v " d 

— (n) =2#(j7SLs,a;(j7SLs))-^(0,a:(0)) + 2 / dr) —${r),x{r))), (115) 

where n is a unit vector in the direction of the observed photons. x(rj) is the 
geodesic along which the light travels from the surface of last scattering (SLS) 
towards us, and ?7sls is the time of last scattering. 

If we choose the topology of the universe to be the orbifold of the Picard 
group, we can use in (112) the Maass cusp forms and the Eisenstein series 
computed in sections 6-9. Let us further choose the initial conditions to be 

f k (0)= - ° ka and f'M = , (116) 

[LBBS97, ITS00, AS01, Lev02], which carry over to a Harrison-Zel'dovich 
spectrum having a spectral index n — 1 and selecting only the non-decaying 
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modes, a is a constant independent of k which is fitted to the amplitude of 
the observed temperature fluctuations. The quantities <j k are random signs, 
a k G {-1,1}. 

The following cosmological parameters are used (2 m = 0.3, Qa = 0.6, J7 C = 
1— J?tot = l-QT-Q m -QA,H{rto) = lOO/iokms^Mpc -1 with h = 0.65. The 
density J? r ~ 10~ 4 is determined by the current temperature To = 2.725 K. 
The point of the observer is chosen to be at x b s = 0.2 + O.li + 1 . 6 j . The 
sky map of figure 12 is computed with the expansion (112) using cusp forms 
and Eisenstein series. For numerical reasons the infinite spectrum is cut such 
that only the eigenvalues with E = k 2 + 1 < 19601 and their corresponding 
eigenfunctions are taken into account. In addition the integration of the con- 
tinuous spectrum is approximated numerically using a Gauss quadrature with 
16 integration points per unit interval. The resulting map in figure 12 is not in 
agreement with the cosmological observations [BHH+03] which are shown in 
figure 14. Especially in the direction of the cusp, the temperature fluctuations 
of our calculated model show a strong peak. Clearly, such a pronounced hot 
spot is unphysical, since it is not observed. This peak comes from the contri- 
bution of the Eisenstein series. The sky map in figure 13, which results from 
taking only the cusp forms into account, yields a much better agreement with 
the cosmological observations. 




Fig. 12. The sky map of the calculated temperature fluctuations of the CMB for 
J?tot = 0.9, On = 0.3, Qa = 0.6, ho = 0.65, and a; obs = 0.2 + O.li + 1.6j, where the 
Eisenstein series and the cusp forms are taken into account. The coordinate system 
is oriented such that the cusp is located in the left part of the figure where the hot 
spot can be seen. 

The hot spot in the direction of the cusp, which results from the contribu- 
tion of the Eisenstein series, is sometimes used, see e.g. [SS76], as an argument 
that the Eisenstein series should not be taken into account. Some papers even 
ignore the existence of the Eisenstein series completely, e.g. [LBBS97]. 
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Fig. 13. The sky map of the calculated temperature fluctuations of the CMB for 
J?tot = 0.9, Q m = 0.3, Q A = 0.6, h = 0.65, and :r obs = 0.2 + O.li + 1.6j, if only the 
cusp forms are taken into account. 




Fig. 14. The temperature fluctuations of the CMB observed by WMAP. Down- 
loaded from [GH03]. 



Another possibility is to take the Eisenstein series into account, but to 
choose for them initial conditions which differ from (116) such that the inten- 
sity towards the cusp is not increasing. To see whether this is possible, let us 
consider the completeness relation of Maass waveforms, 

${r l ,x) = YsWk n ,$M)1>k n {x)+ / ^f scn M V ,-)Hf scn (x)dk, (117) 

n>0 J ° 

cf. (44) and (46), saying that within the fundamental cell any desired metric 
perturbation ^(?7sls, %) a t the time of last scattering can be expressed via the 
conditions 
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fk n (VShs) = W>fc„,#(»7SLS, •)> for thc discrete spectrum, 

fk ' lscn (VSLs) = (i>I lscn , $(VSLS, ')) for the continuous spectrum. 

E.g., choosing the metric perturbation to be 

^(?7SLS,a;) = cos(fc' birrs) (118) 

would neither result in a hot spot in the direction of the cusp nor would 
.ff lscn (v) vanish identically. A more general ansatz would be 

/>oo 

^(vsls,x) = / dk' (A k >(rj SLS ,x) cos(k' In x 2 ) + By (r] S i J s,x)sin(k' In x 2 )), 
Jo 

(119) 

where Ay (?7sls , x ) an d By (77SLS 7 x ) are some amplitudes. 

Concerning the initial conditions we found that (116) has to be modified 
for the continuous part of the spectrum. Hence, the important question occurs: 
How do the correct initial conditions look like? Unfortunately, this question 
has not been answered yet. Work in progress is concerned with finding the 
correct initial conditions and their physical motivation [ALST04]. 

Not knowing the correct initial conditions, we use (116) for the discrete 
part and neglect the continuous part of the spectrum. We proceed with the 
discussion of the properties that come from the discrete part of the spectrum, 
only. 



13 Comparison with the cosmological observations 

The key experiments that observed the fluctuations of the CMB are: COBE, 
Smoot ct al. [SBK+92]; Boomerang, de Bernardis et al. [dBAB+00]; MAXI- 
MA-1, Hanany et al. [HAB+00]; ACBAR, Kuo et al. [KAB+02]; CBI, Pearson 
et al. [PMR+03]; WMAP, Bennett et al. [BHH+03]. 

Concerning the topology of the universe which manifests itself in the low 
multipole moments, there exist only the cosmological observations from COBE 
and WMAP. These two experiments have scanned a full sky map of the tem- 
perature fluctuations. The NASA satellite mission of the Cosmic Background 
Explorer (COBE) was the first experiment that detected the temperature fluc- 
tuations of the CMB in 1992. In addition, it found that the quadrupolc mo- 
ment of the temperature fluctuations is suppressed. In 2003 the NASA satellite 
mission of the Wilkinson Microwave Anisotropy Probe (WMAP) scanned a 
full sky map of the CMB with a much higher resolution confirming and im- 
proving the results of COBE. All the other experiments were either ground 
based or balloon flights that could only scan parts of the sky and thus did 
not measure the first multipole moments. Their strength was to measure the 
finer-scale anisotropy manifested in the higher multipole moments. 
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Fig. 15. The angular power spectrum of the temperature fluctuations observed by 
WMAP. Downloaded from [GH03]. 
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In order to quantitatively compare our results with the observations, we in- 
troduce the angular power spectrum and the correlation function. Expanding 
the temperature fluctuations in the CMB into spherical harmonics, 

oo / 

5T(n(0,<f>)) =J2 E a*m*T(M), (120) 

1=2 m=-l 

yields the expansion coefficients ai m . We remark that in (120) the sum over 
I starts at I — 2, i.e. the monopole and dipole have been subtracted, because 
the monopole term does not give rise to an anisotropy and the dipole term 
is unobservable due to the peculiar velocity of our earth relative to the back- 
ground. From the expansion coefficients we obtain the multipole moments of 
the CMB anisotropics, 

1 ' 

Q ^ 27TT ^ |a ' m|2 ' (121) 

m— — I 

the angular power spectrum, 

and the two-point correlation function, 

CM) = (ST(n)5T(n')) , .„, = —! — -/ f dhdti 5T(n)5T(h'). 

V ' X V ' K '/cos^n-n' 87T 2 sint? J J S 2 xS 2 V ' V ' 

(123) 



n-n —COS 1 



Assuming a Gaussian value distribution for the expansion coefficients aim, the 
two-point correlation function is related to the multipole moments via 

1 oo 

C{^) = — y (21 + 1)0^(008 d). (124) 
47r ^ 

When determining the correlation function numerically, we consider (124) as 
its definition irrespectively whether the expansion coefficients ai, n are Gauss- 
ian distributed or not. 

Figures 15-17 show the angular power spectrum and the correlation func- 
tion, respectively, measured by WMAP [SVP + 03]. In figure 15 we see that 
the first multipole moments resulting from the cosmological observations are 
suppressed, especially the quadrupole, C 2 - The next few multipole moments 
up to I w 25 give rise to a plateau in the angular power spectrum, before the 
larger multipole moments with I > 25 start to increase until they reach the 
first acoustic peak. Page et al. [PNB+03] find that the first acoustic peak is 
at I = 220.1 ± 0.8. The trough following this peak is at I = 411.7 ± 3.5, and 
the second peak is at I = 546 ± 10. 
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Fig. 16. The angular power spectrum of the calculated temperature fluctuations 
for l?tot = 0.9, Sl m = 0.3, Q A = 0.6, h = 0.65, and x obs = 0.2 + O.li + 1.6j (□) in 
comparison with the angular power spectrum of the WM AP observations ( x ) . 



Figures 16 and 17 show the angular power spectrum and the correlation 
function, respectively, corresponding to the calculated sky map for the Pi- 
card model of figure 13 in comparison with the results of the cosmological 
observations and with the concordance model [SVP + 03]. Regarding our re- 
sults of the calculated temperature fluctuations, see figure 16, we see that the 
first multipole moments are suppressed in accordance with the cosmological 
observations. Especially the quadrupole term, C2, computed in our model, 
is suppressed. However the suppression of the observed one is still stronger. 
Note, that in the case of the concordance model, the first multipoles are in- 
creased in contrast to our results. This points to a non-trivial topology of our 
universe. 

The next multipole moments give rise to a plateau in the angular power 
spectrum. We use this plateau to fit the constant a of the initial conditions 
(116) such that the sum of the first 19 computed multipole moments matches 
with the cosmological observations, 

20 20 

^ ' Cj, computed — ^ ' C(, observed- (1^5) 
Z=2 Z=2 
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Fig. 17. The solid line is the correlation function of the calculated temperature 
fluctuations for l? to t = 0.9, Q m = 0.3, Qa = 0.6, ho = 0.65, and :r b s = 0.2 + O.li + 
1.6j. The dashed line represents the correlation function of the WMAP observations, 
and the dotted line shows the concordance model [SVP + 03]. 



But concerning the fluctuations on finer scales (smaller angular resolution), 
i.e. higher multipoles, we do not obtain the acoustic peaks in our numerical 
calculations. This is due to the fact that we neglect several physical effects 
which perturb the temperature on finer scales. We ignore all the physical 
effects that dominate on scales smaller than those corresponding to I ^> 25, 
since our main interest is in the influence of the non-trivial topology of the 
universe due to the Picard group. 

Another comparison of our calculated temperature fluctuations with the 
cosmological observations can be done via the correlation function (124) which 
emphasises large scales, i.e. the behaviour of low multipoles. Concerning the 
correlation function C(?9), we find quite good agreement of our calculated 
temperature fluctuations with the cosmological observations, see figure 17. 
The reason for this good agreement of our model with the observations is 
due to the suppressed quadrupole moment. We also show in figure 17 (dotted 
line) the concordance model [SVP + 03] which is not in good agreement with 
the data for > 7° in contrast to our model (full curve). Especially for 
large angular separations, $ > 160°, the concordance model is not able to 
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describe the observed anticorrelation in C($). This anticorrelation constitutes 
a fingerprint in the CMB that favours a non-trivial topology for the universe. 

14 Conclusion 

We pointed out the importance of the eigenvalue equation of the Laplacian 
and gave some physical applications. In the framework of quantum chaos, we 
focused on the eigenvalue equation of the Picard orbifold and computed the 
solutions numerically. Our main goal was to determine the statistical proper- 
ties of the solutions, but also to apply the solutions to a completely different 
problem, namely the temperature fluctuations in the CMB. Here, we demon- 
strate that a model of the universe with the topology of the Picard group 
matches the large-scale anisotropy in the CMB better than the current con- 
cordance model. Thereby, we show that the methods developed in quantum 
chaos can be successfully applied to other fields of research although the phys- 
ical interpretation can differ completely. 
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